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A solution of codimension 2 brane is found for which 4 dimensional Friedmann cosmology is 
recovered on the brane with time dependent tension, in the Einstein frame. The effective parameter 
p/p of equation of state on the brane can be quintessence like, de Sitter like or phantom like, 
depending on integration constants of the solution. 
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In 1983 Rubakov and Shaposhnikov presented an extraordinary model pj to solve the cosmological constant problem, 
I/"") ' in which the cosmological constant appeared as an integration constant. In this model the effective cosmological 
constant was free of the zero point quantum energy. The zero point energy only warps the extra dimension, and 
at the same time our 4 large dimensional spacetime keeps flat. This is similar to the case that although in realistic 
universe the 4 dimensional curvature may be high but the 3 dimensional spatial curvature is almost zero. 

During recent years codimension 1 brane model has been fairly well understood. Also the cosmology in frame of 
codimension 1 brane has been thoroughly studied, for a review see Q. Generally speaking codimension 1 cosmology 
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of RS model implies Friedmann equation like 



H 2 ~p+ P -, (1) 
a 

> 

where a is the brane tension, p is the energy density of matter on the brane. We find the standard cosmology recovers 
when the energy scale is low enough and the corrected Friedmann equation gives notable effect when the energy scale 
is high enough . But to the codimension 2 case, we get into a much more weird situation: although the effective 
cosmological constant is independent of the tension of the brane, to our surprise, we appear to be very restricted in 
1^ 1 our permitted brane energy-momentum. Typically, a brane in its ground state has a very special energy-momentum 
tensor, which is isotropic and has the property that energy density=-pressure Q, that is to say, the only spacetime 
allowed is maximally symmetric space, no Friedmann cosmology with generic parameter of state equation p/p can be 
set up on the brane. 

To circumvent this embarrassed situation and get a viable cosmology on codimension 2 brane, different authors have 
put much effort. Clearly, two possible ways are available. First, inspired by the early work on cosmic string y|, a fat 
brane is a natural way to circumvent this problem. A such model has been considered in in which a codimension 
2 braneworld with spherical extra dimensions compactified by magnetic flux is studied. Assuming Einstein gravity, 
it shows that when the brane contains matter with an arbitrary parameter of equation of state, and find that the 
universe expands consistently with standard Friedmann cosmology. However, the cost is the relation between the 
brane tension and the bulk deficit angle becomes 5 ~ p — 3p for a general equation of state. This relation does 
not imply a self-tuning of the effective 4D cosmological constant. So this progress is bankrupt in the merit that the 
effective cosmological constant is independent of the brane tension, which is not what we want. In Q the action 
includes classical Hilbert action, we may curious whether or not a higher order Ricci scalar terms can circumvent the 
restriction energy dcnsity=-prcssurc on the brane. That is the second accessible possibility. A proposal along this 
direction has been presented in 0, in which a general braneworld in 6 dimensional Einstein-Gauss-Bonnet gravity 
is studied. It shows how the 4 dimensional Einstein equation is recovered for the induced metric and matter on the 
brane. It also shows that relaxing regularity of the curvature in the vicinity of the brane gives rise to an additional 
possible correction to the Einstein equation. Now within this proposal the 4 dimensional Einstein equation recovers 
on the brane, therefore the Friedmann cosmology can be obtained on the brane without question. But we may 
bring forward the following questions to this suggestion: why in pure classical domain we must consider higher order 
correction to revive Einstein equation on the brane? If the R 2 term correction gives the answer that we need, as the 
authors stressed, since 0(R 2 ) corrections to the Hilbert Lagrangian do arise in the low energy limit of string theory, 
the inclusion of this type of term could be regarded as mandatory if one wants to embed any braneworld solution into 
string/M-theory, how about the more higher order corrections? They may give "false" results. 
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In this paper we successfully recover almost realistic Friedmann cosmology on the brane in frame 6 of dimensional 
Einstein gravity with a 3 brane on which the tension is variable. If the brane tension originates from the quantum 
zero energy, it would decay with the expansion of the universe, as suggested in 8] . That article addresses the question 
of whether non-perturbative effects of self-interacting quantum fields in curved space-times may yield a significant 
contribution to the tension. Focusing on the trace anomaly of quantum chromo-dynamics (QCD), a preliminary 
estimate of the expected order of magnitude yields 

p = nH, (2) 

where n = m 3 , m is energy scale of QCD chiral phase transition, p is vacuum energy (brane tension) and H is Hubble 
parameter. H is defined as 

H = - d i- w 

a at 

As a natural generalization of Zi symmetry in codimension 1 brane model, the metric ansatz is taken to be 

ds 2 = -dt 2 + a{t) 2 dx 2 + dr 2 + b{t) 2 v(r) 2 de 2 , (4) 

where dx 2 is the metric of 3 dimensional maximally symmetric space, v(r) = sin(r), r, or, sinh(r), depending on 
the internal space is a sphere, plane or pseudosphere. To spherical case the internal space is compact by nature. And 
to case of plane or pseudosphere the internal space is compactfied by proper identification, that is to say, we identify 
all the points r =constant. The brane stands at r = 0. Supposing the metric is regular on the brane, it is easy to get 
the familiar result in codimension 2 infinite thin braneworld models 0, Q 

1 - b = 4G 6 cr, (5) 

where a is the brane tension which equals the energy density of the brane, and Gq denotes the 6 dimensional Newton 
constant. According to Schutzhold's investigation the brane tension is proportional to the Hubble parameter, as 
shown in equation J3J). In its original form the coefficient is related to the energy scale of QCD chiral phase transition. 
Although it gives right energy scale of our observed cosmological constant, we should not restrict us to the energy 
scale of QCD phase transition from the theoretical view because various phase transitions may happen in the history 
of our universe. Here we adopt a more generic case, 

p(t) = a(t) = -±-cH, (6) 

in which c is a constant, it can contain all the contributions of phase transitions ever happened on the brane. Note 
that the dimension of c is [I]. Thus a smaller c corresponding to a higher energy scale. Einstein field equations of this 
system is 

^G AB = 8ttG 6 T ab , (7) 

where 

T AB = T^+T^ B \ (8) 

^Gab denotes 6 dimensional Einstein tensor, Tab is the total energy momentum tensor, and T A b g represent 

the energy on the brane and in the bulk respectively, and Gq is the 6 dimensional Newton constant. The brane energy 
momentum is shown clearly in matrix form, 

fS - ( f! *' Slr> o | . (9) 

where T^J^ = —cr(t)g^ v , and p or v runs from to 3. We do not try to specify here because we will find later 

that it has to take a form of non perfect fluid for a self consistent solution under condition ©. Now consider the 
functions a(t) and b(t) in the following forms 



2 ( c 2 +t) 2 (£ + c 2 )V4 + 2c 2 ci 
a(t) = c 3 e " cosh [ — ], (10) 
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and 



(t + c 2 )V4 + 2c 2 ci. 



m = i vr^tohi ^^^ 01 ]. (ii) 

Here ci, C2, C3 are 3 integration constants. We shall explain their physical meanings later in Einstein frame. One can 
verify that a(t) and b(t) solve the 6 dimensional Einstein equation with non perfect fluid in the bulk. We can prove 
it is the unique solution satisfying the following conditions: 

1. The metric is spatially flat, that is, dx 2 in is Euclidean. 

2. Although we permit the matter in the bulk to be non perfect fluid, we impose the condition — T^ bu ' = T^- bu \ — 
rp(bu)^ _ £ Q ensure t na t it behaves as perfect fluid along the brane. From the physical viewpoint it is a 
reasonable requirement. 

Here we will present some other notes on this solution. Both a(t) and b(t) are internal-curvature-independent, that is, 
they do not depend on the Gaussian curvature of the internal space, which means, spherical, flat and pseudospherical 
internal spaces correspond to the same a(t) in l|10(l . The tension of the brane is time dependent, which most probably 
leads to the time dependence of the energy momentum tensor in the bulk. We find it is just the case. Explicitly, 

n 3cosh[/]- 3 ((4-cic 2 )cosh[/] + (12 + 5cic 2 )cosh[/] + 2^(4+ (4+ ^ coshf/1) sinhf/])) 
_ T (bu) v _ V f / / 12 \ 

16c 2 7rG 6 (l + Uanh[/]) ' V ' 

and 

T (bu)* = T (6„)i + _i_ f^2 cosh[/] -2 + 6(4 + CiC 2 + 2/ tanh[/])^) , (13) 

where / = l Ss^l and I = ^/A + 2 Cl c 2 . -T (bu *>° = T (bu *>\ = T (bu )\ = T^l and are also internal-curvature- 

independent. The situation is different to T*- 6 "^, it is internal-curvature-dependent. To avoid complexity we only 
give specific form of the flat internal space case, 

, 72 + 30cic 2 - ^cosh[/]" 2 + 6^ + U + c lC 2 - cos\i[f]~ 2 \ tanh[/] 

rp(bu) _ rp(bu) j f V f / 

5 1 87rG 6 c 2 (l + Ztanh[/]) ' V ' 

To the cases of sphere and pseudosphere the discussions are similar. This model also possesses another remarkable 
property. There is a non zero energy flux along the r direction emerging in the bulk. In the case of flat internal space, 
it has the form as 

,„ _ rp{bu) _ 12 cosh[/]- 2 

W4 ~ 1 40 ~ 167rG 6 cr(l + Ztanh[/])- (15) 

This is exactly a rational result. We see the vacuum energy of the brane is decaying, but there is nothing on the 
brane to absorb the energy decayed |9| . Therefore the only way out of this energy is to flow into the bulk. There 
is a singularity in (|15fl at r = 0, where the brane habitats. Any finite energy flux flowing from the brane must be 
divergent at the neighborhood of the brane , since the brane is infinite thin. In respect to these results we are able to 
predict the energy flux vanishes at r = tt in case of spherical internal space and slopes more sharply with respect to 
r in case of pseudospherical internal space. In fact these are just the main differences of bulk energy fluxes between 
flat and spherical or pseudospherical internal spaces. Still another point deserves to address. There is a deficit angle 
in the extra dimension due to the existence of the brane in all previous works on the codimension 2 brane models 
0010] Here, however, from fTTJl it is clear the "deficit angle" is negative, or, there is a surplus angle 27r(|6| — 1) 
since |6| > 1. For example the internal space becomes a mandarin orange rather than an American football in case of 
spherical internal space. 

If we believe only Einstein frame makes sense in our real world, it is not obvious whether or not the solution 1|10|) 
really represents a generic Friedmann universe. Now we turn to Einstein frame, in which the metric becomes 

ds% = ~\b\dt 2 + \b\a 2 dx 2 . (16) 

Define proper time 



dr = \b\^ 2 dt, 



(17) 
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and the proper scale factor, 

d=\b\^ 2 a. (18) 

Define cosmic velocity 

da 

= & (19) 



and cosmic acceleration 



A.-g. (20, 



Now it is time for explaining the physical meanings of the constants c\, C2, C3. Obviously C3 is only a constant scale 
factor, and C2 represents a relocation of initial time in (|1U|) . To make out the meaning of ci, we calculate the proper 
Hubble parameter 

H=^. (21) 

We find 

H(t = -c 2 ) = - + ^. (22) 
c 2 

It is clear that in some sense c\ indicates initial cosmic velocity. Naively when c 2 c\ = —2 this solution degenerates 
to a de Sitter brane without angle deficit in the bulk in view of b = — 1. Even this most simple case is not the trivial 
one studied before, because in our case a is not a constant all the same, while a =constant in all present studies. 

The most significant parameters from the viewpoint of observations are parameter of state equation w — cnC rgy S dcnsity 
and the deceleration parameter q of the total cosmic fluid. The two parameters are closely related to each other. 
Therefore, we calculate the effective w and q corresponding to the proper scale factor and proper time. By definition 

PcS V 2 + ^ A c / O o% 
Weff = = ^75 , (23) 

Pes 2>V£ 

and 



aA c 

v e 

we arrive at 

4cic 2 (3 + cic 2 ) cosh[/] - (10 + § I 2 + Z 4 ) cosh[/] + 3(4 + 3c lC 2 + 2 Cl cH sinh[/] - 2(3 + Cl c 2 )l sinh[/]) 
q= (ci - 2(3 + cic 2 ) cosh[/] - 3/ sinh[/]) 2 u ' 

and 



(24) 



(25) 



Wctf = ~- + -q, (26) 



where u = y/l + I tanh[/]. The equations are rather complicate. It needs some labour to check when c 2 ci > —2, 
the universe finally becomes quintessence like phase, which in some sense can simulate dark matter plus cosmological 
constant universe; when c 2 c\ — —2, it is a de Sitter phase; and when c 2 c\ < —2, the universe goes into a phantom like 
phase. This is our main result: Friedmann universe of rather arbitrary parameters of state equations can recover on 
the codimension 2 varying-tension-brane. Different types of the universe depend on different choices of the integration 
constants. We show this result clearly in some figures. The figures illustrate that the evolution of the universe only 
depends on c 2 c\. Recall equation 10, c is the vacuum energy on the brane. It is clear we can always choose a proper 
ci to counter the effect of an arbitrarily large c. Accordingly the most important virtue of prior models with 2 extra 
dimensions 0, 0| is retained in our model. One may be afraid that if c 2 ci < —2, the arguments of the hyperbolic 
functions in l|25|l become imaginary numbers. But we find the final results always keep real even if the imaginary 
arguments appear. 
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FIG. 1: - versus w e s and q: the left part. 

To showing clearly the properties of w e s and q we plot two figures with the same set of parameters, in which t is in different 
intervals of real line. To this figure - 6 (0,0.3), and to the next | £ (0.3,3). In this figure cic 2 = 1 , 02 has been tuned such 
that w e ff = 1 when t — 0. 
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FIG. 2: I versus w e s and q: the right part. 

In this figure the parameters are chosen as same as the last figure. These two figures show that a universe undergoes stiff 
matter, radiation, dust like phases and finally becomes a quintessence like universe. From the observations we know that w e g 
is declining, so we may think that the current universe corresponds to the coordinate - = 0.3. If it is true in the future w will 
decrease for some time and then increase, eventually reach approximately -0.68. 

Because the age of our universe is model dependent we can not obtain zero energy from | = 0.3 by using the value of t in the 
standard model. But we point out that for any given c we can always obtain such a universe, as shown in figure I and figure 
2, by tuning ci . 





FIG. 3: - versus w. 

c 

In this figure Cic 2 = —0.5. This figure shows that we can adopt parameters to ensure the final value of w equals 0.72, which 
may offer some implications on the coincidence problem. 
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FIG. 4: | versus u>. In this figure cic 2 



2.1. A phantom dominated universe emerges under this parametrization. 



To sum up: in this paper we discuss a codimcnsion 2 braneworld model, in which the brane tension depends on time. 
We find a unique solution under some reasonable considerations. According to this solution Fricdmann cosmologies , 
which permit different parameters of state equations, can be set up on the brane in Einstein frame. This is sharply 
different from the previous researches on this topic. At the same time our model keeps the most impressive character of 
the codimension 2 brane model — the effective cosmological constant is free of the absolute value of the brane tension. 
In our model the parameter of state equation only depends on the integration constants, that is to say, it circumvents 
the famous 10 120 contradiction. So the problem becomes why the integration constants are taking those values. 

We will also point out some questions to be discussed future in our model. Follow from figure ^ if it is believed 
that we are living at t/c — 0.3, there exists another coincidence problem, which is essentially the choices of integration 
parameters, all the same. The asymptotic behavior of the universe in figure |3] overcomes the coincidence problem, but 
it is difficult to simulate the history of our universe, especially the history of structure formation. It is still an open 
problem in present stage. 
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